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In the extraction of the charm contribution F£ c to the proton structure function F 2 in our re- 
cent publication (T), we have not properly taken into account the running of the electromagnetic 
coupling a cm . The measured cross sections were corrected to the Born level for QED radiation, 
but not for the running of a ern . This was not taken properly into account in the extraction of 

Z7CC 
b 2 ■ 

In addition, the cross section predictions of the CASCADE program were calculated with 
fixed a em . The cross section in the visible range calculated with running a em is 5.63 nb (instead 
of 5.09 nb given in (TJ). The conclusions on the description of the data by CASCADE are 
unchanged. The extrapolation factors, defined as the ratio of the full cross section cj^ to 
the cross section cr^ s co in the visible phase space of the D* meson, and their uncertainties are 
changed slightly. The amended values are shown in figure [TJ which replaces figure 15 of (T). 

The amended values of F 2 CC extracted from measured D* ± cross sections with the HVQDIS 
program and with the CASCADE program are lower by about 6 upto 11% as compared to [QJ. 
The corrected values of F 2 CC and its uncertainties are given in table [TJ which replaces table 1 1 
of (TJ. The amended F 2 CC values are compared to a measurement based on lifetime information 
determined with the HI silicon vertex detector [2] and with theoretical predictions in figures [2] [3] 
and[4l which replace figures 16, 17 and 18 of [1], respectively. 
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Table 1: F 2 CC in bins of Q 2 and x extracted from measured D* cross sections with two dif- 
ferent programs, HVQDIS and CASCADE. The extrapolation uncertainty 5 ext is determined 
by varying model parameters within a program. The statistical (5 sta t) and systematic (S sys t) 
uncertainties arise from the determination of the D* cross section and are the same for both 
programs. 



2 



o 

£ .2 

V 

0) — 

£ 3 



[_ Q 2 = 6.5 GeV 2 j 

- rn 

z 

' I ' 


: Q 2 = 12.0 GeV 2 

I 


[_ Q 2 = 20.0 GeV 2 j 


: Q 2 = 35.0 GeV 2 

: I 
I 

■ ' 


Q 2 = 60.0 GeV 2 j 

" i i i . ." 


10" 4 10" 3 10* 

X 

a CASCADE (AO) 
• HVQDIS(CT10f3) 



10" 4 10" 3 10" 2 



Figure 1 : Extrapolation factors from the visible phase space to the total phase space for the D* 
meson as determined from HVQDIS and CASCADE. The error bars show the extrapolation 
uncertainty which is determined by varying the theory parameters listed in tables 1 and 2 of [HI. 
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Figure 2: F 2 CC as derived from D* data with HVQDIS (points). The inner error bars show the 
statistical uncertainty, the outer error bar the statistical and experimental systematic uncertainty 
added in quadrature. The extrapolation uncertainty within the HVQDIS model is shown as 
blue band in the bottom of the plots. The outer (orange) band shows the model uncertainty 
obtained from the difference in F 2 CC determined with HVQDIS and CASCADE. The data are 
compared to the measurement of F 2 CC with the HI vertex detector [2] (open squares), to NLO 
DGLAP predictions from HVQDIS with two different proton PDFs, and to the F| c prediction 
of HERAPDF1.0. 
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Figure 3: F 2 CC as derived from D* data with HVQDIS (points). The inner error bars show the 
statistical uncertainty, the outer error bar the statistical and experimental systematic uncertainty 
added in quadrature. The extrapolation uncertainty within the HVQDIS model is shown as blue 
band in the bottom of the plots. The outer (orange) band shows the model uncertainty obtained 
from the difference in F| c determined with HVQDIS and CASCADE. The data are compared 
to the measurement of F| c with the HI vertex detector [2] (open squares) and to predictions 
from the global PDF fits CT10 (dashed line), MSTW08 at NNLO (dark dashed-dotted line), 
NNPDF2.1 (shaded band) and ABKM (light dashed-dotted line). 
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Figure 4: F| c as a function of Q 2 for different x, as derived from D* data with HVQDIS (points). 
The inner error bars show the statistical uncertainty, the outer error bar the total uncertainty, 
including statistical, experimental systematic, extrapolation and model uncertainty added in 
quadrature. The data are compared to the measurement of F 2 CC with the HI vertex detector [|2) 
(open squares), to NLO DGLAP predictions from HVQDIS with two different proton PDFs, 
and to the Ff prediction of HERAPDF1 .0. 
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Abstract 

Inclusive production of D* mesons in deep-inelastic ep scattering at HERA is stud- 
ied in the range 5 < Q 2 < 100 GeV 2 of the photon virtuality and 0.02 < y < 0.7 of 
the inelasticity of the scattering process. The observed phase space for the D* meson is 
Pt{D*) > 1.25 GeV and \r](D*)\ < 1.8. The data sample corresponds to an integrated 
luminosity of 348 pb _1 collected with the HI detector. Single and double differential cross 
sections are measured and the charm contribution F| c to the proton structure function F2 is 
determined. The results are compared to perturbative QCD predictions at next-to-leading 
order implementing different schemes for the charm mass treatment and with Monte Carlo 
models based on leading order matrix elements with partem showers. 



Accepted by Eur. Phys. J. C 



F. D. Aaron 5,48 , C. Alexa 5 , V. Andreev 25 , S. Backovic 30 , A. Baghdasaryan 38 , 

S. Baghdasaryan 38 , E. Barrelet 29 , W. Bartel 11 , K. Begzsuren 35 , A. Belousov 25 , P. Belov 11 , 
J.C. Bizot 27 , M.-O. Boenig 8 , V. Boudry 28 , 1. Bozovic-Jelisavcic 2 , J. Bracinik 3 , G. Brandt 11 , 
M. Brinkmann 11 , V. Brisson 27 , D. Britzger 11 , D. Bruncko 16 , A. Bunyatyan 13 ' 38 , 

G. Buschhorn 26 ^, L. Bystritskaya 24 , A.J. Campbell 11 , K.B. Cantun Avila 22 , F. Ceccopieri 4 , 
K. Cerny 32 , V. Cerny 16 ' 47 , V. Chekelian 26 , J.G. Contreras 22 , J.A. Coughlan 6 , J. Cvach 31 , 
J.B. Dainton 18 , K. Daum 37 ' 43 , B. Delcourt 27 , J. Delvax 4 , E.A. De Wolf 4 , C. Diaconu 21 , 

M. Dobre 12 ' 50 ' 51 , V. Dodonov 13 , A. Dossanov 26 , A. Dubak 30 - 46 , G. Eckerlin 11 , S. Egli 36 , 
A. Eliseev 25 , E. Elsen 11 , L. Favart 4 , A. Fedotov 24 , R. Feist 11 , J. Feltesse 10 , J. Ferencei 16 , 

D. -J. Fischer 11 , M. Fleischer 11 , A. Fomenko 25 , E. Gabathuler 18 , J. Gayler 11 , S. Ghazaryan 11 , 

A. Glazov 11 , L. Goerlich 7 , N. Gogitidze 25 , M. Gouzevitch 11 - 45 , C. Grab 40 , A. Grebenyuk 11 , 
T. Greenshaw 18 , B.R. Grell 11 , G. Grindhammer 26 , S. Habib 11 , D. Haidt 11 , C. Helebrant 11 , 
R.C.W. Henderson 17 , E. Hennekemper 15 , H. Henschel 39 , M. Herbst 15 , G. Herrera 23 , 

M. Hildebrandt 36 , K.H. Hiller 39 , D. Hoffmann 21 , R. Horisberger 36 , T. Hreus 4 ' 44 , F. Huber 14 , 
M. Jacquet 27 , X. Janssen 4 , L. Jonsson 20 , A.W. Jung 15 , H. Jung 11,4 ' 52 , M. Kapichine 9 , 
I.R. Kenyon 3 , C. Kiesling 26 , M. Klein 18 , C. Kleinwort 11 , T. Kluge 18 , R. Kogler 11 , P. Kostka 39 , 
M. Kraemer 11 , J. Kretzschmar 18 , K. Kriiger 15 , M.P.J. Landon 19 , W. Lange 39 , 

G. Lastovicka-Medin 30 , P. Laycock 18 , A. Lebedev 25 , V. Lendermann 15 , S. Levonian 11 , 
K. Lipka 11,50 , B. List 12 , J. List 11 , R. Lopez-Fernandez 23 , V. Lubimov 24 , A. Makankine 9 , 

E. Malinovski 25 , P. Marage 4 , H.-U. Martyn 1 , S.J. Maxfield 18 , A. Mehta 18 , A.B. Meyer 11 , 

H. Meyer 37 , J. Meyer 11 , S. Mikocki 7 , 1. Milcewicz-Mika 7 , F. Moreau 28 , A. Morozov 9 , 
J.V. Morris 6 , M. Mudrinic 2 , K. Miiller 41 , Th. Naumann 39 , PR. Newman 3 , C. Niebuhr 11 , 

D. Nikitin 9 , G. Nowak 7 , K. Nowak 11 , J.E. Olsson 11 , D. Ozerov 24 , P. Pahl 11 , V. Palichik 9 , 

I. PanagouhW' 11,42 , M. Pandurovic 2 , Th. Papadopoulou' 1 11,42 , C. Pascaud 27 , G.D. Patel 18 , 

E. Perez 10 - 45 , A. Petrukhin 11 , 1. Picuric 30 , S. Piec 11 , H. Pirumov 14 , D. Pitzl 11 , R. Placakyte 12 , 

B. Pokorny 32 , R. Polifka 32 , B. Povh 13 , V. Radescu 14 , N. Raicevic 30 , T. Ravdandorj 35 , 
P. Reimer 31 , E. Rizvi 19 , P. Robmann 41 , R. Roosen 4 , A. Rostovtsev 24 , M. Rotaru 5 , 

J.E. Ruiz Tabasco 22 , S. Rusakov 25 , D. Salek 32 , D.P.C. Sankey 6 , M. Sauter 14 , E. Sauvan 21 , 
S. Schmitt 11 , L. Schoeffel 10 , A. Schoning 14 , H.-C. Schultz-Coulon 15 , F. Sefkow 11 , 
L.N. Shtarkov 25 , S. Shushkevich 26 , T. Sloan 17 , 1. Smiljanic 2 , Y. Soloviev 25 , P. Sopicki 7 , 
D. South 11 , V. Spaskov 9 , A. Specka 28 , Z. Staykova 11 , M. Steder 11 , B. Stella 33 , G. Stoicea 5 , 
U. Straumann 41 , T. Sykora 4 ' 32 , PD. Thompson 3 , T. Toll 11 , TH. Tran 27 , D. Traynor 19 , 
P. Truol 41 , 1. Tsakov 34 , B. Tseepeldorj 35 - 49 , J. Tumau 7 , K. Urban 15 , A. Valkarova 32 , 

C. Vallee 21 , P. Van Mechelen 4 , Y. Vazdik 25 , D. Wegener 8 , E. Wiinsch 11 , J. Zacek 32 , 
J. Zalesak 31 , Z. Zhang 27 , A. Zhokin 24 , H. Zohrabyan 38 , and F. Zomer 27 

1 /. Physikalisches Institut der RWTH, Aachen, Germany 

2 Vinca Institute of Nuclear Sciences, University of Belgrade, 1 100 Belgrade, Serbia 

3 School of Physics and Astronomy, University of Birmingham, Birmingham, Uff 

4 Inter-University Institute for High Energies ULB-VUB, Brussels and Universiteit Antwerpen, 
Antwerpen, Belgium 

5 National Institute for Physics and Nuclear Engineering (NIPNE) , Bucharest, Romania" 1 

6 Rutherford Appleton Laboratory, Chilton, Didcot, UK b 

7 Institute for Nuclear Physics, Cracow, Poland d 

8 Institut fiir Physik, TU Dortmund, Dortmund, Germany 11 

9 Joint Institute for Nuclear Research, Dubna, Russia 

10 CEA, DSM/Irfu, CE-Saclay, Gif-sur-Yvette, France 

11 DESY, Hamburg, Germany 



1 



12 InstitutfUr Experimentalphysik, Universitat Hamburg, Hamburg, Germany 11 

13 Max-Planck-Institut ftir Kernphysik, Heidelberg, Germany 

14 Physikalisches Institut, Universitat Heidelberg, Heidelberg, Germany 11 

15 Kirchhoff- Institut fur Physik, Universitat Heidelberg, Heidelberg, Germany a 

16 Institute of Experimental Physics, Slovak Academy of Sciences, Kosice, Slovak Republic? 

17 Department of Physics, University of Lancaster, Lancaster, UK b 

18 Department of Physics, University of Liverpool, Liverpool, UK b 

19 Queen Mary and Westfield College, London, UK b 

20 Physics Department, University of Lund, Lund, Sweden 9 

21 CPPM, Aix-Marseille Univ., CNRS/IN2P3, 13288 Marseille, France 

22 Departamento de Fisica Aplicada, CINVESTAV, Merida, Yucatan, Mexico^ 

23 Departamento de Fisica, CINVESTAV IPN, Mexico City, Mexico^ 

24 Institute for Theoretical and Experimental Physics, Moscow, Russia k 

25 Lebedev Physical Institute, Moscow, Russia 6 

26 Max-Planck-Institut fur Physik, Munchen, Germany 

27 LAL, Universite Paris-Sud, CNRS/IN2P3, Orsay, France 

28 LLR, Ecole Polytechnique, CNRS/IN2P3, Palaiseau, France 

29 LPNHE, Universite Pierre et Marie Curie Paris 6, Universite Denis Diderot Paris 7, 
CNRS/IN2P3, Paris, France 

30 Faculty of Science, University of Montenegro, Podgorica, Montenegro 11 

31 Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republic 11 

32 Faculty of Mathematics and Physics, Charles University, Praha, Czech Republic 11 

33 Dipartimento di Fisica Universita di Roma Tre and INFN Roma 3, Roma, Italy 

34 Institute for Nuclear Research and Nuclear Energy, Sofia, Bulgaria 6 

35 Institute of Physics and Technology of the Mongolian Academy of Sciences, Ulaanbaatar, 
Mongolia 

36 Paul Scherrer Institut, Villigen, Switzerland 

37 Fachbereich C, Universitat Wuppertal, Wuppertal, Germany 

38 Yerevan Physics Institute, Yerevan, Armenia 

39 DESY, Zeuthen, Germany 

40 InstitutfUr Teilchenphysik, ETH, Zurich, Switzerland 1 

41 Physik-Institut der Universitat Zurich, Zurich, Switzerland 1 

42 Also at Physics Department, National Technical University, Zografou Campus, GR-15773 
Athens, Greece 

43 Also at Rechenzentrum, Universitat Wuppertal, Wuppertal, Germany 

44 Also at University of P.J. Safdrik, Kosice, Slovak Republic 

45 Also at CERN, Geneva, Switzerland 

46 Also at Max-Planck-Institut fiir Physik, Munchen, Germany 

47 Also at Comenius University, Bratislava, Slovak Republic 

48 Also at Faculty of Physics, University of Bucharest, Bucharest, Romania 

49 Also at Ulaanbaatar University, Ulaanbaatar, Mongolia 

50 Supported by the Initiative and Networking Fund of the Helmholtz Association (HGF) under 
the contract VH-NG-401. 

51 Absent on leave from NIPNE-HH, Bucharest, Romania 

52 On leave of absence at CERN, Geneva, Switzerland 
t Deceased 

a Supported by the Bundesministerium fUr Bildung und Forschung, FRG, under contract 



2 



numbers 05H09GUF, 05H09VHC, 05H09VHF, 05H16PEA 

b Supported by the UK Science and Technology Facilities Council, and formerly by the UK 
Particle Physics and Astronomy Research Council 

c Supported by FNRS-FWO-Vlaanderen, IISN-IIKW and IWT and by Interuniversity Attraction 
Poles Programme, Belgian Science Policy 

d Partially Supported by Polish Ministry of Science and Higher Education, grant 
DPN/N168/DESY/2009 

e Supported by the Deutsche Forschungsgemeinschaft 

f Supported by VEGA SR grant no. 2/7062/27 

9 Supported by the Swedish Natural Science Research Council 

h Supported by the Ministry of Education of the Czech Republic under the projects LC527, 

INGO-LA09042 and MSM002 1620859 

1 Supported by the Swiss National Science Foundation 

j Supported by CONACYT, Mexico, grant 48778-F 

k Russian Foundation for Basic Research (RFBR), grant no 1329.2008.2 and Rosatom 

1 This project is co-funded by the European Social Fund (75%) and National Resources (25%) 

- (EPEAEK II) - PYTHAGORAS II 

m Supported by the Romanian National Authority for Scientific Research under the contract 
PN 09370101 

n Partially Supported by Ministry of Science of Montenegro, no. 05-1/3-3352 



3 



1 Introduction 



The measurement of the charm production cross section and the derived structure function F 2 CC 
in deep-inelastic electrorfl-proton scattering (DIS) at HERA allows tests of the theory of the 
strong interaction, quantum chromodynamics (QCD). Previous measurements ifD-HTl of charm 
production in DIS at HERA have demonstrated that charm quarks are predominantly produced 
by the boson gluon fusion process, which is sensitive to the gluon density in the proton. The 
production of charm quarks contributes up to 30% to the inclusive ep scattering cross section. 
The correct treatment of effects related to the charm quark contribution in perturbative QCD 
calculations, in particular the mass effects, is therefore important for the determination of parton 
distribution functions (PDFs). 

At HERA several different techniques have been used to determine the charm contribu- 
tion F| c to the proton structure function F 2 . Besides the full reconstruction of a D or D* 
meson JMIIMSIIMIl, the lifetime of heavy flavoured hadrons [T7ti9l[T2l[T4l[T7l or the semi- 
leptonic decay [fT3l are exploited. Compared to the other methods, the measurement of D* 
mesons provides a charm sample with a large signal-to-background ratio. The results presented 
here are based on a data sample collected by the HI experiment, corresponding to an integrated 
luminosity of 348 pb -1 . Increased statistics, extended phase space, as well as reduced system- 
atic uncertainties compared to the previous HI analysis [[TOl make more detailed tests of pQCD 
predictions possible. Compared to earlier HI analyses the phase space for the D* meson is 
extended in transverse momentum from p T (D*) > 1.5 GeV to p T (D*) > 1.25 GeV and in 
pseudo-rapidity from \r](D*)\ < 1.5 to \i](D*)\ < 1.8. This extension reduces the amount of 
extrapolation needed for the determination of F| c . 



2 QCD Models and Monte Carlo Simulation 

The QCD models employed for data corrections and for comparison with measured cross sec- 
tions are introduced in the following. Different Monte Carlo (MC) generators based on leading 
order (LO) QCD calculations complemented with parton showers are used to simulate detector 
effects in order to determine the acceptance and the efficiency for selecting DIS events with a 
D* meson and to estimate the systematic uncertainties associated with the measurements. The 
generated events are passed through a detailed simulation of the HI detector response based on 
the GEANT program [18] and through the same reconstruction and analysis algorithms as used 
for the data. 

The MC program RAPGAP [19] is used for the generation of the direct process of boson 
gluon fusion to a cc pair. It uses a LO matrix element with massive charm quarks. Parton show- 
ers [|20l based on the DGLAP evolution equations 11211 model higher order QCD effects. The 
hadronisation of partons is performed with PYTHIA [|22l which implements the Lund String 
Model [|23l . For the fragmentation of the charm quark into the D* meson the Bowler param- 
eterisation [|2~4l is chosen and the longitudinal part of the fragmentation function is reweighted 
to the Kartvelishvili parameterisation ll25Tl . The parameter a of the Kartvelishvili parameter- 
isation is set to the values measured by HI [26] which depend on the centre-of-mass energy 
squared s of the hard subprocess (yg — > cc). The threshold between the two regions in s 

'in this paper "electron" is used to denote both electron and positron. 
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is chosen such that the mean value of s in the lower region is in agreement with the mean s 
of the event sample without a jet associated with the D* meson [|26l . RAPGAP is interfaced 
to the HERACLES program [27| in order to simulate the radiation of a real photon from the 
incoming or outgoing lepton and virtual electro-weak effects. For the determination of the de- 
tector acceptance and efficiency, RAPGAP is used with the PDF set CTEQ6.6M J28l which 
is derived at next-to-leading order (NLO), but gives a good description of the data. Alterna- 
tively, RAPGAP is used with CTEQ6LL ||29ll derived at LO. The mass of the charm quark is 
set to m c = 1.5 GeV. The renormalisation scale [i r and the factorisation scale /ij are set to 
[i r = jif = \J Q 2 + 4m;? + (p T ) 2 , where Q 2 denotes the photon virtuality and p T the trans- 
verse momentum of the charm quark in the photon-gluon centre-of-mass frame. The relevant 
parameter settings and their variations are summarised in tabled] 



RAPGAP 


Parameter name 


Central value 


Variation 


Charm mass 


m c = 1.5 GeV 




Renormalisation scale 


Mr = VQ 2 + 4m c + (Pt) 2 




Factorisation scale 


n = VQ 2 + Am2 c + (PT) 2 




Fragmentation 


a = 10.3 for s < s threshold 

a = 4.4 for S > Sthreshold 
Sthreshold = 70 GeV 


8.7 < a < 12.2 
3.9 < a < 5.0 

50 < sthreshold < 90 GeV 2 


PDF 


CTEQ6.6M 


CTEQ6LL 


CASCADE 


Parameter name 


Central value 


Variation 


Charm mass 


m c = 1.5 GeV 




Renormalisation scale 


AV,o = VQ 2 + 4m c + Pt 


1/2 < HrjUrfS < 2 


Factorisation scale 


A*/,o = \f s + Qt 


1/2 < Hf/nffl < 2 


Fragmentation 


a = 8.4 for S < S threshold 
a = 4.5 for S > S threshold 
Sthreshold = 70 GeV 


7.3 < a < 9.7 
3.9 < a < 5.1 

50 < threshold < 90 GeV 2 


PDF 


AO 


jjt r variation: AO-, A0+ 



Table 1: Parameters used in the MC simulations. The central choice of the renormalisation 
(factorisation) scale is denoted by /i r (/i/,o)- The invariant mass squared and the transverse 
momentum squared of the cc pair are denoted by s and Q\, respectively, m c is the charm quark 
mass and p* T and p T are the transverse momentum of the charm quark in the photon-gluon 
centre-of-mass frame and in the electron-proton centre-of-mass frame, respectively, a is the 
fragmentation parameter in the Kartvelishvili parameterisation. Two values of a in two regions 
of s with the boundary ^threshold are used [1261 . For CASCADE different PDF sets are available 
which were determined for a variation of the renormalisation scale by a factor of 1/2 or 2. These 
are used consistently for the // r variation here. 
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The CASCADE [ 30] program is based on A; T -factorisation and the CCFM evolution equa- 
tions [31 1. In CASCADE the direct boson gluon fusion process is implemented using off-shell 
matrix elements convolved with a /c r -unintegrated gluon distribution of the proton. The PDF set 
AO [1321 is used. Time-like parton showers of the charm quark and anti-quark are implemented, 
but those from initial state gluons are not. The hadronisation of partons is performed in the same 
way as for RAPGAP. When CASCADE is used for the extrapolation to F 2 CC , the renormalisation 
and factorisation scales are varied to estimate the theoretical uncertainty. For the variation of 
the renormalisation scale, the PDF sets AO- and A0+ are used, which were extracted with the 
corresponding scale variation [32] . The parameter variations used in CASCADE are also listed 
in table [Q 

In addition to RAPGAP and CASCADE, the data are also compared to predictions of an 
NLO calculation |[33l[34ll based on collinear factorisation and the DGLAP evolution equa- 
tions. This calculation assumes three active flavours (u, d, s) in the proton (fixed-flavour-number 
scheme: FFNS) and massive charm quarks are produced dynamically via boson gluon fusion. 
The predictions are calculated using the program HVQDIS [|34l . Corresponding fixed-flavour 
NLO parton density functions of the proton, CT10f3 [1331 (with the strong coupling set to 
a s (M z ) = 0.106) and the NLO variant of MSTW2008f3 (351, are used. Charm quarks are 
fragmented in the ^p centre-of-mass frame into D* mesons using the Kartvelishvili |[25l pa- 
rameterisation for the fragmentation function with the value of the parameter a as measured by 
HI 11261 . The renormalisation and factorisation scales are set to ji r — pif — v/Q 2 + 4m 2 .. The 
value used for the charm mass is 1.5 GeV. To obtain the theoretical systematic uncertainty for 
the extrapolation to F 2 CC the parameters are varied according to table |2] Each of the variations 
is performed independently. The resulting uncertainties are added in quadrature to obtain a 
conservative estimate of the total theoretical uncertainty. 

The results are also compared with a NLO calculation [38] based on the zero-mass variable- 
flavour-number scheme (ZM-VFNS), where the charm quark is considered as a massless con- 
stituent of the proton. This calculation is only valid for a sufficiently large transverse momentum 
of the D* meson p^{D*) in the 'jp centre-of-mass frame. For the comparison to this prediction 
the analysis is therefore restricted to p^(D*) > 2 GeV. The ZM-VFNS uses the fragmentation 
function determined in [|3~9l . The scales are chosen to be // r = /// = a/ (Q 2 + (p^) 2 )/2. 

3 HI Detector 

A detailed description of the HI detector can be found elsewhere 14011411 . Only the components 
essential to the present analysis are described here. The origin of the HI coordinate system is 
the nominal ep interaction point. The direction of the proton beam defines the positive 2-axis 
(forward direction). Transverse momenta are measured in the x-y plane. Polar (9) and az- 
imuthal (0) angles are measured with respect to this reference system. The pseudo-rapidity is 
defined as r] = — mtan(0/2). 

Charged particles are measured within the central tracking detector (CTD) in the pseudo- 
rapidity range —1.85 < rj < 1.85. The CTD consists of two large cylindrical jet chambers 
(CJCs), surrounding the silicon vertex detector CST [|42l . The CJCs are separated by a drift 
chamber which improves the z coordinate reconstruction. A multiwire proportional cham- 
ber [1431 . which is mainly used in the trigger, is situated inside the inner CJC. These detec- 
tors are arranged concentrically around the interaction region in a magnetic field of 1.16 T. 
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HVQDIS 


Parameter name 


Central value 


Variation 


Charm mass 


77i c = 1.5 GeV 


1.3 <m c < 1.7 GeV 


Renormalisation scale 


Mr o = \/Q 2 + 


1/2 < Hr/llrfl < 2 


Factorisation scale 


u f n = JQ 2 + Am 2 


1/2 < Uflutn < 2 


Fragmentation 


a = 6.1 for S < S threshold 
a = 3.3 for S > S threshold 
Sthreshold = 70 GeV 


5.3 < a < 7.0 
2.9 < a < 3.7 
50 < hhreshoid < 90 GeV 2 


PDF 


CT10O 


MSTW2008f3 


Fragmentation fraction 


/(c D*) = 23.8 ± 0.8% El 




ZM-VFNS 


Parameter name 


Central value 


Variation 


Charm mass 


m c = 1.5 GeV 




Renormalisation scale 


A*r,o=vW + (P50 2 )/ 2 


1/2 < Hr/Hrfl < 2 


Factorisation scale 


/*/,o = \/(Q 2 + (Pt) 2 )/2 


1/2 < /i///i/, < 2 


Fragmentation 


KKKS08 f39H 




PDF 


CTEQ6.6M 





Table 2: Parameter variations used for the uncertainty estimation of the NLO calculations. The 
central choice of the renormalisation (factorisation) scale is denoted by //^o (A*/,o)- m c is the 
charm quark mass and a is the fragmentation parameter in the Kartvelishvili parameterisation. 
In the two regions of s, separated by the boundary s t hreshoid, different values of a are used fl26l . 

The trajectories of charged particles are measured with a transverse momentum resolution of 
c(Pt)/pt ~ 0.5%pT/GeV © 1.5% [|44|. The interaction vertex is reconstructed from CTD 
tracks. The CTD also provides triggering information based on track segments measured in 
the CJCs Il45ll46l and a measurement of the specific ionisation energy loss dE/dx of charged 
particles. The forward tracking detector measures tracks of charged particles at smaller polar 
angle (1.5 < 77 < 2.8) than the central tracker. 

Charged and neutral particles are measured in the liquid argon (LAr) calorimeter, which 
surrounds the tracking chambers and covers the range —1.5 < 77 < 3.4 with full azimuthal 
acceptance. Electromagnetic shower energies are measured with a precision of a(E)/E = 
12%/y/E/GeV © 1% and hadronic energies with a{E)/E = 50%/ \J E /GeV © 2%, as de- 
termined in test beam measurements [47] . A lead-scintillating fibre calorimeter (SpaCal) pT| 
covering the backward region —4.0 < 77 < —1.4 completes the measurement of charged and 
neutral particles. In this analysis the SpaCal is used in particular for the identification and re- 
construction of the scattered electron. For electrons a relative energy resolution of a(E)/E = 
7%/ ■sjE I GeV © 1% is achieved, as determined in test beam measurements [48 1 . The SpaCal 
provides energy and time-of-flight information used for triggering purposes. A Backward Pro- 
portional Chamber (BPC) in front of the SpaCal is used to improve the angular measurement of 
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the scattered electron. 

The hadronic final state is reconstructed using an energy flow algorithm which combines 
charged particles measured in the CTD and the forward tracking detector with information 
from the SpaCal and LAr calorimeters [[491. 

The luminosity determination is based on the measurement of the Bethe-Heitler process 
(ep — > ep'j) where the photon is detected in a calorimeter located at z = —103 m downstream 
of the interaction region in the electron beam direction. 



4 Event Selection and Signal Extraction 

The data sample corresponds to an integrated luminosity £ = 348 pb _1 and was recorded with 
the HI detector in e + p (185 pb -1 ) and e~p interactions (163 pb _1 ) in the years 2004 to 2007. 
During this period electrons at an energy of 27.6 GeV were collided with protons at 920 GeV. 
The events were triggered by a local energy deposit in the SpaCal in coincidence with at least 
one track in the CTD, with an overall trigger efficiency of 98%. 

DIS events are selected by requiring a high energy electromagnetic cluster in the SpaCal 
which is consistent with resulting from the scattered electron. The event kinematics including 
the photon virtuality Q 2 , the Bjorken scaling variable x and the inelasticity variable y are re- 
constructed with the eS method ||50~1 , which uses information from the scattered electron and 
the hadronic final state and provides good resolution in the covered y range. The kinematic 
region for the photon virtuality is restricted to 5 < Q 2 < 100 GeV 2 corresponding to the geo- 
metric acceptance of the SpaCal. In order to ensure a high trigger efficiency, the energy of the 
electron candidate is required to fulfil E' e > 10 GeV. The inelasticity is restricted to the range 
0.02 < y < 0.7. 

Charm production is identified by the reconstruction of D* mesons using the decay channel 
D *± _^ jjo^ _^ K^^-nf which has a branching fraction BIZ = 2.63 ± 0.04% [EQ. The 
tracks of the decay particles are reconstructed in the CTD. The invariant mass of the _ft' =F 7r ± 
system is required to be consistent with the D° mass fl5TTl within ±80 MeV. A loose particle 
identification criterion is applied to the kaon candidates using the measurement of the specific 
energy loss, dE/dx, in the CTD. This improves the signal-to-background ratio, especially at 
low transverse momenta of the D* meson. The kinematic range of the measurement is sum- 
marised in table [3] Details of the selection are described in 11521 . 

D*^ candidates are selected using the mass difference method [1531 . The distribution of the 
mass difference Am = m{K T 'K ± 'Kf) — m{K T Ti ± ) is shown in figure [TJ A clear signal peak 
around the nominal mass difference of 145.4 MeV [|5T1 is observed. 

The wrong charge combinations, defined as K ±r K ± ixf with K ± n ± pairs in the accepted 
D° mass range, are used to constrain the shape of the combinatorial background in the signal 
region. The number of D* mesons is determined in each analysis bin by a simultaneous fit 
to the right and the wrong charge Am distribution. As the signal has a tail towards larger 
Am values, the asymmetric Crystal Ball function 11541 is used for the signal description. The 
shape of the background is parameterised with the Granet function [|55l . An unbinned extended 
log likelihood fit Il56l is performed using the RooFit framework 11571 . 
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The fit to the complete data set yields 24705±343 D* mesons. This represents an increase in 
statistics of an order of magnitude compared to the previous analysis IfTOl . For each analysis bin 
the fit to the Am spectrum uses the two parameters describing the signal asymmetry obtained 
from the fit to the complete data set. The width of the peak varies in dependence of the D* 
kinematics and is therefore treated as a free parameter of the fit. 

Photon virtuality Q 2 I 5 < Q 2 < 100 GeV 2 

Inelasticity y 0.02 < y < 0.7 

Pseudo-rapidity of D m± -1.8 < 7](D*) < 1.8 

Transverse momentum of D*^ 1 pt(D*) > 1.25 GeV 

Table 3: Definition of the kinematic range of the present measurement. 



5 Cross Section Determination and Systematic Errors 

The following formula is used to calculate the inclusive D* meson production cross section at 
the Born level in the visible kinematic range defined in table [3] 

a vis (ep -> eD*X) = ^^^TTXa S " (1) 

£ • BK{D* -> Kti-k s ) • (1 + <5 rad ) 

Here N(D* ± ) is the number of D* mesons obtained using an unfolding procedure defined 
below, r is the contribution from reflections from other decay modes of the D° meson, C is the 
integrated luminosity, BIZ is the branching ratio and 5 rad denotes the QED radiative corrections. 
For the differential measurements the cross section is also divided by the bin width. No bin 
centre correction is applied. 

To obtain N(D*) in each measurement bin, the data are corrected for detector effects in- 
cluding the trigger efficiency by means of regularised matrix unfolding I158H6H . The response 
matrix A which reflects the acceptance and the resolution of the HI detector relates the distri- 
butions y rcc of reconstructed variables to distributions x true of variables at the generator level, 
Ax true = y vec . Each matrix element Ay is the probability for an event originating from bin 
j of x tTue to be measured in bin % of y rec . The response matrix is determined from simulation 
and has twice as many bins at the reconstruction level as at the generator level in order to pro- 
vide detailed information on the probability distribution and to improve thereby the accuracy 
of the unfolding procedure. The procedure reduces statistical correlations between neighbour- 
ing bins and the influence of model assumptions in the cross section determination. Additional 
bins outside of the kinematic range of this measurement are used to provide constraints on the 
migrations into the measured phase space. 

The measured D* cross section includes decays from B hadrons to D* mesons which are 
expected to contribute to less than 2%. For the determination of F 2 CC the beauty contribution as 
calculated with HVQDIS is subtracted. 

For the present analysis the systematic error dominates over the statistical uncertainty for 
almost the whole phase space. The measurement is statistically limited only for large transverse 
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momenta or large photon virtualities. The systematic uncertainties are determined in each bin 
separately and are summarised in table |4] for the total cross section. They are divided into 
uncertainties which are bin-to-bin uncorrected and uncertainties which are correlated between 
the bins. The uncertainties in the following are given in percent of the cross section values. 

The following sources for bin-to-bin uncorrected systematic errors are considered: 

Signal Extraction: Using different parameterisations for the signal and background shapes Il52~| 
the systematic uncertainty on the signal extraction is estimated to be 2%. 

Radiative Corrections: For the correction of the measured cross section to the Born level, the 
HERACLES interface to RAPGAP is used. The corrections are of the order of 2.5%. An 
uncertainty of 2% is assigned ifTTl . 

Trigger efficiency: The efficiency of the trigger conditions requiring an energy deposition in 
the SpaCal and a central track is at least 95%. The combined uncertainty is estimated to 
be 1%. 

D° mass cut: The invariant mass resolution of the data is not fully reproduced by the MC sim- 
ulation, leading to different efficiencies of the D° mass cut. The difference is evaluated 
by comparing the width of the D° signal in data and MC and extrapolating to the region 
outside of the mass cut assuming a Gaussian distribution [52J. The dependence on the D* 
kinematics is studied, and the maximum difference of 1.5% is assigned as uncertainty. 

Reflections: The amount of reflections r from decay modes of the D° meson other than D° — > 
K~^7i^~ is determined using a RAPGAP MC sample of inclusive charm events and is 
found to amount to 3.8%. The dependence of r on kinematic quantities is studied in the 
simulation and found to be constant within 1%, a value which is used as the systematic 
uncertainty. 

Photoproduction background: The photoproduction background is estimated using a PYTHIA 
photoproduction MC sample to be less than 0.2%, which is used as systematic uncertainty. 

dE/dx cut: The loss of D* signal events due to the dE/dx requirement on the kaon track 
amounts to 3.4% in data. The dE/dx cut is not applied in the simulation, but corrected 
for in the data by a global factor. The dependence of the cut efficiency on kinematic 
variables is studied and found to be within 2%, which is used as systematic uncertainty. 

Where appropriate, the effect of the bin-to-bin correlated systematic uncertainties is eval- 
uated by changing the response matrix and repeating the unfolding procedure. The following 
correlated error sources are considered: 

Track finding efficiency: The systematic error on the track efficiency of 4.1% per D* meson 
is the dominant error of this measurement. It arises from two contributions: (i) The 
comparison of the track finding efficiency in data and simulation leads to an error of 2% 
for the slow pion track and 1% for the tracks of the D° daughter particles and is assumed 
to be correlated between the decay particles; (ii) the efficiency with which a track can be 
fitted to the event vertex leads to a systematic error of 1% per D* meson. The uncertainty 
on the track finding efficiency is considered to be half correlated between the bins of the 
measurement. 
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Luminosity: The systematic error on the luminosity measurement is estimated to be 3.2%. 

Branching ratio: The uncertainty on the branching ratio of the D* meson is 1.5% [j5TTl . 

Model: The parton shower model uncertainty is on average 2%, estimated by taking the differ- 
ence in cross section obtained using RAPGAP or CASCADE for the data correction. 

PDF: Using different parton density functions in RAPGAP for the data correction leads to an 
uncertainty of below 1%. 

Fragmentation: The parameter of the Kartvelishvili fragmentation function in RAPGAP is 
varied in the range given in table [TJ The resulting differences in the cross section are 
between 1% and 5%. 

Electron energy: The systematic uncertainty on the SpaCal energy scale is 0.5% which results 
in a systematic error of typically below 1%, but up to 10% at large D* inelasticity z (see 
section©. 

Electron angle: The angular resolution of the SpaCal/BPC of 0.5 mrad leads to a systematic 
uncertainty of typically 2%. 

Hadronic energy: The systematic uncertainty on the energy scale of the hadronic final state is 
2%. The influence in general is small (below 0.5%) but leads to larger uncertainties of up 
to 20% at large D* inelasticity z(D*) and small y. 

All sources of the systematic errors are assumed to be uncorrected between the sources and 
added in quadrature. This results in an overall systematic uncertainty of 7.6%. 



6 Cross Sections 

The cross section in the visible range defined in table |3]is measured to be: 

a vis (ep eD^X) = 6.44 ± 0.09 (stat.) ± 0.49 (syst.) nb . (2) 

The corresponding predictions from RAPGAP (CTEQ6LL), RAPGAP (CTEQ6.6M), and CAS- 
CADE (AO) amount to 5.02 nb, 4.37 nb, and 5.09 nb, respectively. The NLO calculation 
HVQDIS yields 5.98+^ nb with CT10f3 as the proton PDF and 5.52+°;^ nb with MSTW2008f3, 
where the uncertainty is determined by varying the parameters according to tableland adding 
the resulting uncertainties in quadrature. The HVQDIS predictions are slightly below the mea- 
surement but agree with the data within errors. 

In tableland figures |2] and |3] the single differential cross sections are presented as a func- 
tion of variables describing the event kinematics: the photon virtuality Q 2 , the inelasticity y 
and Bjorken x. The measurements are compared to the predictions of the MC programs RAP- 
GAP and CASCADE (figure |2]) and of the next-to-leading order calculation with the HVQDIS 
program (figure 13} . Since the theoretical calculations predict smaller cross sections than the 
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Uncorrelated uncertainties 




& /o 


Radiative correction 


2% 


Trigger efficiency 


1% 


D° meson mass cut 


1.5% 


Reflections 


1.0% 


Photoproduction background 


< 0.2% 


dE/dx cut 


2% 


Correlated uncertainties 


Track efficiency 


4.1% 


Luminosity 


3.2% 


Branching ratio 


1.5% 


Model 


2.1% 


PDF 


1% 


Fragmentation 


2.6% 


Electron energy scale 


1.3% 


Electron angle 9 


1.3% 


Hadronic energy scale 


0.3% 


Total systematic uncertainty 


7.6% 



Table 4: Summary of all sources of systematic uncertainties and their effect on the D* pro- 
duction cross section with the breakdown into bin-to-bin uncorrelated and bin-to-bin correlated 
sources. 

measurement, the normalised ratio _R norm of theory to data is shown in the lower part of the fig- 
ures, which facilitates the shape comparison between the different theoretical predictions and 
the data. This ratio is defined as: 

i / theo . do- tho ° 

morrn / vis ' dY /i\ 

~~ i /-data da data ' ^ ' 

x / u vis dY 

where and dcr theo /dF are the total and differential cross section of the model under con- 
sideration and Y denotes any measured variable. The normalisation uncertainties of the data 
(luminosity, branching ratio and half of the tracking uncertainty) cancel in this ratio. Similarly, 
uncertainty sources of the NLO predictions altering only the normalisation do not affect R noim 
since for each variation the total and the differential cross section are varied simultaneously. 
In all predictions the decrease with Q 2 is slightly less steep than in data. The y dependence is 
reasonably well described by all predictions. The dependence on x is slightly less steep in all 
predictions than in the data, an effect which is larger for the NLO calculations. 

In tableland figures |4]and|5]the single differential cross sections are presented as a function 
of the kinematic variables of the D* meson: the transverse momentum pt(D*) and pseudo- 
rapidity r)(D*) in the laboratory frame, the transverse momentum p^(D*) in the 7p centre- 
of-mass frame and the D* inelasticity z(D*). The D* inelasticity z(D*) corresponds to the 
fraction of the virtual photon momentum carried by the D* meson and is determined as z(D*) = 
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(P ■ Pd*)I (P ■ q) = (E — p z ) d* /2yE e where E e is the energy of the incoming electron and P, 
q and p D * denote the four-momenta of the incoming proton, the exchanged photon and the D* 
meson, respectively. All predictions are able to describe the shape of the p T (D*) distribution of 
the data reasonably well, although RAPGAP has a tendency to underestimate the data at large 
Pt(D*). The shape of the r] distribution shows sensitivity to the use of different parton densities 
in the RAPGAP MC. The prediction based on CTEQ6.6M agrees better with the data than the 
prediction based on CTEQ6LL. A very good description of the r] shape is obtained with the 
CASCADE MC. The HVQDIS calculations with CT10f3 and MSTW2008f3 both describe the 
7] distribution reasonably well, but have a tendency to be too low in the positive 77 (forward) 
region. For the transverse momentum p^(D*) in the ^p centre-of-mass frame, the RAPGAP 
MC using either PDF is too steep at large p* T {D*\ while the CASCADE prediction generally 
has a different shape. The NLO predictions are in good agreement with the data. The z(D*) 
dependence is not described by any of the calculations, showing a deficit of all predictions at 
low z(D*) values. 

In order to investigate the correlation between pseudo-rapidity and transverse momentum, a 
double differential measurement in p T (D*) and t](D*) is performed. In tableland figures [6]and 
|7]the double differential cross section is presented as a function of the pseudo-rapidity i](D*) in 
bins of the transverse momentum of the D* meson pt(D*). In the backward direction almost 
no D* mesons with large transverse momentum are produced. At low transverse momenta 
all predictions are below the data in the very forward direction. At pt(D*) > 6 GeV the 
CASCADE and HVQDIS predictions give a good description of the data, while RAPGAP is 
too low. 

While the transverse momentum of the D* meson in the laboratory frame is correlated with 
the photon virtuality Q 2 , the transverse momentum in the ^p centre-of-mass frame p* T {D*) is 
directly related to the hard subprocess. The double differential cross section as a function of 
r](D*) and p^(D*) are presented in tableland figures [8] and [9] The distribution shows simi- 
lar behaviour to the double differential cross section as a function of r](D*) and the transverse 
momentum p T (D*) in the laboratory frame. They are in general better described by the predic- 
tions of CASCADE and HVQDIS, while RAPGAP underestimates the data for positive rj(D*) 
at large p* T (D*). 

The double differential cross section measurements in y and Q 2 are presented in table [8] 
and figures \T0\ and \TT\ All predictions are able to describe the distribution reasonably well, 
independent of the PDF choice. At low Q 2 RAPGAP as well as CASCADE has a tendency to 
be lower than the data. 

To allow a comparison to the ZM-VFNS predictions, the cross sections are also measured 
with an additional cut p^(D*) > 2 GeV. The measurements are shown in tables [5] and [6] and in 
figures Hi] and [13] The ZM-VFNS calculation overshoots the data at low y. The x dependence 
is less steep than for the data, and it has a different shape in pr(D*). The dependence of the 
cross section on the other variables is described reasonably well. In general the ZM-VFNS 
prediction describes the data worse than the NLO FFNS calculation HVQDIS. Also at higher 
Q 2 > 100 GeV 2 the ZM-VFNS prediction fails to describe the D* production lfT5l . 

In order to facilitate the comparison with previous measurements the cross sections are also 
measured in a reduced phase space of the D* meson: p T (D*) > 1.5 GeV and \r](D*)\ < 1.5. 
They are listed in tables [9] and [TOj In figure [14] these measurements are shown as a function 
of Q 2 together with the results of the measurement at high Q 2 [fT5l . These measurements span 
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over almost three orders of magnitude in Q 2 . The data are well described by CASCADE and 
the HVQDIS predictions with both PDF sets in the whole Q 2 range, while RAPGAP overshoots 
the data at high Q 2 . 



7 Charm Contribution to the Proton Structure Function 

The charm contribution F£°(x, Q 2 ) to the proton structure function F 2 (x, Q 2 ) is related to the 
charm cross section in the one photon exchange approximation by: 

<1V< 2 ^ : ' n ([l + (l-yf]Ff(x,Q 2 )-y 2 Fl%x,Q 2 )) . (4) 



dxdQ 2 Q 4 x 

Weak interaction effects are neglected and a em denotes the electromagnetic coupling constant. 
The contribution from the structure function F[ c is less than 4% in the present x — Q 2 range. 
Assuming the ratio F£ C /F 2 CC is predicted correctly within a model, the measured inclusive 
cross sections a^(y, Q 2 ) in bins of y and Q 2 are converted to a bin centre corrected 
F| cexp ((a;), (Q 2 )) using the relation Q 2 = xys and extrapolating cr vis to the full phase space: 



Ff exp ((x), (Q 2 » = °£ (y ' Q *[ ■ Ff theo ((x), (Q 2 )) . (5) 

[y,Q 



Here and F£ ctheo are the theoretical predictions. The HVQDIS program is used to cal- 
culate cr^ s co and F 2 CC theo in the NLO DGLAP scheme. In the kinematic range of the current 
analysis the beauty contribution to the D* cross section is small. It is estimated with HVQDIS 
and subtracted for the determination of F 2 CC . 

The measurement covers about 50% of the total phase space for charm production. The 
extrapolation factor to the full phase space is model dependent. Since CASCADE also provides 
a reasonable description of the cross sections in the phase space covered in this analysis, it is 
used as an alternative model to determine F 2 CC exp in order to investigate this model dependence. 
The extrapolation factors in the present analysis, defined as the ratio of the full cross section 
to the cross section in the visible phase space of the D* meson, determined with 
HVQDIS and CASCADE, are compared in figure [15] They differ by about 10% at medium x. 
In the largest x bin the extrapolation factor as well as the difference between the two models 
increases significantly. In general the extrapolation factor determined with HVQDIS has smaller 
uncertainties than the one from CASCADE. Due to the larger phase space of the D* meson 
in the present analysis compared to previous measurements, where the phase space coverage 
amounted to about 30%, the extrapolation factor to the full phase space is considerably smaller, 
although the model dependence of the extrapolation remains sizable. 

In table [TT] and figure [16] the resulting F 2 C extracted from the inclusive cross sections 
with HVQDIS is shown as a function of x for different values of Q 2 . In addition to the exper- 
imental systematic uncertainties described in section |6] the extrapolation (equation |5]) leads to 
an uncertainty. This extrapolation uncertainty is determined by varying the theory parameters 
listed in table [2] simultaneously in the calculation of cr^ s eo and F 2 CC theo . The resulting uncertain- 
ties on F 2 C are shown separately in figure [161 HVQDIS and CASCADE both give a reasonable 
description of the measured cross sections and can be used to extract F 2 CC . The differences ob- 
tained in F 2 CC exp for the two models are used to define the model uncertainty on F 2 CC , which is 
also given in figure [ 
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The results of a F| c measurement based on lifetime information determined with the HI 
silicon vertex detector CST IfPfl is compared to the present measurement in figure [161 The two 
measurements are based on independent methods with similar precision and agree very well. 
The F 2 CC determined with D* mesons covers a larger range in x due to the larger r\ coverage of 
the CJCs compared to the CST It also has smaller uncertainties at low Q 2 , where the uncertainty 
of the lifetime based measurement is dominated by the light quark background. 

Figure [16] also compares the FFNS NLO calculation of F 2 c to the measurement using the 
MSTW2008B and CT10f3 proton PDFs. Both calculations give a reasonable description of 
the data. The measurement is also compared to the F 2 CC prediction for HERAPDF1.0 II621 , 
which has been extracted from the HI and ZEUS combined inclusive proton structure function 
data. A general-mass variable-flavour-number scheme [|63l has been used which interpolates 
between the FFNS and the ZM-VFNS. The uncertainty on Ff for the HERAPDF1.0 predic- 
tion is dominated by the variation of the charm mass in the PDF fit, which is included in the 
model uncertainty of the prediction. In general the prediction agrees with the F 2 CC measurement, 
showing that the gluon density determined from the scaling violations of the inclusive DIS cross 
section is consistent with that observed in charm production. At large x the central value of the 
HERAPDF1 .0 prediction has a tendency to lie above the F 2 C data, which may indicate a prefer- 
ence for a larger charm mass than the central value used for HERAPDF1.0, m c = 1.4 GeV. 

The F 2 CC measurement is compared to predictions from global PDF fits in figure[T7] CT10 [|35l , 
MSTW2008 NNLO (64l and NNPDF2.1 [65J have been derived in general-mass variable- 
flavour-number schemes, while for the ABKM fit ll66l the FFNS including higher order ra- 
diative corrections in QCD adopting the running MS mass has been used. All predictions give 
a reasonable description of the measurement. At low Q 2 all predictions have a tendency to 
decrease less steeply with x than the data. 

The measured F 2 CC as a function of Q 2 for different values of x is shown in figure [T_8] Scaling 
violations are observed. The Q 2 dependence of the data is well reproduced by the FFNS NLO 
calculation, but at low x the predictions are below the data, an effect which is larger if the 
MSTW2008O set is used as proton PDF. The HERAPDF1.0 prediction is in agreement with 
the data. 



8 Conclusion 

A measurement of D* meson production in deep-inelastic scattering is performed with a tenfold 
increase in statistics and a significantly enlarged phase space compared to the previous HI 
measurement. Single and double differential cross sections are determined as a function of 
variables describing the kinematics of the event as well as of the D* meson. The measurements 
are found to be reasonably well described by predictions based on the fixed-flavour-number 
scheme, namely the leading order Monte Carlo simulations RAPGAP and CASCADE as well 
as the next-to-leading order calculation HVQDIS. The data are also compared to a next-to- 
leading order calculation in the zero-mass variable-flavour-number scheme, which in general 
describes the data less well and is particularly high at low y. 

The double differential cross section as a function of Q 2 and y is used to determine the 
charm contribution F 2 C to the proton structure function F 2 . The extrapolation to the full phase 
space is done with two different models, using the next-to-leading order calculation HVQDIS 
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and the Monte Carlo program CASCADE based on leading order matrix elements and parton 
showers. The results for F 2 CC in these two models are very similar except for the highest x 
values. The results agree well with a measurement based on lifetime information determined 
with the HI vertex detector. The data are well described by next-to -leading order calculations 
using different PDFs, showing that the gluon density determined from the scaling violations of 
the inclusive DIS cross sections is consistent with the one observed in charm production. 
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Q 2 range [GcV 2 ] 


dcr/dQ 2 [nb/GeV 2 ] 


da/dQ 2 [nb/GeV 2 ] 
tor p f T (D*) > 2.0 GeV 


5.0 
6.0 


6.0 
8.0 


0.782 
0.538 


±0.048 
±0.022 


+0.058 
0.057 
+0.039 
=0.039 


0.317 
0.237 


±0.023 
±0.012 


+0.023 
-=0.022 
+0.018 
==0.018 


8.0 


10.0 


0.384 


±0.018 


+0.028 
=0.028 


0.1568 


±0.0094 


+0.0115 
==0.0110 


10.0 


13.0 


0.249 


±0.011 


+0.018 

=o.oi8 


0.1156 


±0.0063 


+0.0085 
==0.0084 


13.0 


19.0 


0.1549 


±0.0057 


+0.0109 

=0.0110 


0.0695 


±0.0031 


+0.0051 
=0.0050 


19.0 


27.5 


0.0874 


±0.0038 


+0.0062 
=0.0062 


0.0350 


±0.0019 


+0.0025 
=0.0025 


27.5 


40.0 


0.0463 


±0.0022 


+0.0032 
=0.0032 


0.0206 


±0.0011 


+0.0015 
=0.0015 


40.0 
60.0 


60.0 
100.0 


0.0188 
0.00824 


±0.0013 
±0.00057 


+0.0014 
=0.0014 
+0. 00058 
=0.00057 


0.00856 
0.00478 


±0.00067 
±0.00037 


+0.00059 
=0.00059 
+0.00034 
=0.00034 


y range 


dcr/dy [nb] 


dcr/dy [nb] 
lor p^(D*) > 2.0 GeV 


0.02 


0.05 


21.67 


±1.04 


+2.53 
=2.83 


4.75 


±0.40 


+0.81 
=0.65 


0.05 


0.09 


20.97 


±0.94 


+1.51 
= 1.79 


7.14 


±0.41 


+0.68 
=0.71 


0.09 


0.13 


20.05 


±0.97 


+ 1.48 

=i.58 


7.61 


±0.46 


+0.68 
=0.65 


0.13 
0.18 


0.18 
0.26 


14.63 
12.61 


±0.80 
±0.54 


+1.04 
= 1.03 
+0.91 

=0.90 


6.79 
5.01 


±0.50 
±0.25 


+0.52 
=0.52 
+0.39 
=0.38 


0.26 


0.36 


8.39 


±0.43 


+0.72 
=0.63 


4.25 


±0.20 


+0.32 
=0.32 


0.36 


0.50 


5.87 


±0.31 


+0.63 
=0.47 


2.96 


±0.17 


+0.25 
=0.24 


0.50 


0.70 


3.00 


±0.27 


+0.32 
=0.30 


1.83 


±0.16 


+0.19 
=0.18 


x range 


dcr/dx [nb] 


dcr/dx [nb] 
forp5,(D*) > 2.0 GeV 


0.00007 


0.00020 


4990 


±300 


+440 

^go 


2970 


±180 


+230 
=210 


0.00020 


0.00035 


6020 


±280 


_i_460 
=c 440 


3060 


±160 


+220 
=210 


0.00035 


0.00060 


4180 


±170 


_i_320 
^310 


1994 


±93 


+ 143 
= 141 


0.00060 


0.00100 


2631 


±109 


_i_190 

^iss 


1172 


±57 


+87 
=86 


0.00100 


0.00170 


1540 


±61 


_i_107 

^ios 


586 


±31 


+42 
=42 


0.00170 


0.00330 


579 


±24 


±tf 


235 


±12 


+22 
=22 


0.00330 


0.05000 


13.24 


±0.61 


_i_1.17 
^1.33 


4.20 


±0.24 


+0.45 
=0.45 



Table 5: Differential D* cross section as a function of Q 2 , y and x in the kinematic range of 
5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7, \rj(D*)\ < 1.8 and p T (D*) > 1.25 GeV. The first 
quoted uncertainty is statistical and the second is systematic. 
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Pt range [GeV] 



dcr/dp T [nb/GeV] 



da/dp T [nb/GeV] 
forp* (£>*) > 2.0 GeV 



1.25 


1.60 


2.55 


±0.21 


±0 


1.60 


1.88 


2.88 


±0.19 


i 
±0 


1.88 


2.28 


2.68 


±0.11 


i o 


z.zo 


Z.DO 


z.±4( 


±U.UoD 


4-0 


2.68 


3.08 


1.538 


±0.058 


=t 


3.08 


3.50 


1.362 


±0.047 


±° 


3.50 


4.00 


0.924 


±0.032 


=t 


4.00 


4.75 


0.630 


±0.020 




4.75 


6.00 


0.2987 


±0.0098 


±° 


6.00 


8.00 


0.0883 


±0.0039 


=t 


8.00 


11.00 


0.0217 


±0.0015 


±° 


11.00 


20.00 


0.00183 


±0.00034 


=t 



0.334 


±0.042 


±° 


0.436 


±0.055 




0.853 


±0.057 




0.935 


±0.057 


4- 


0.744 


±0.047 


^O 


0.806 


±0.042 


4- 


0.620 


±0.033 


=t 


0.443 


±0.022 


± 


0.239 


±0.012 


4- 


0.0769 


±0.0042 


=t 


0.0210 


±0.0016 


=t 


0.00188 


±0.00032 


± 



.035 
.035 
.041 
042 
.061 
.062 
.086 
.086 
.065 
.065 
.061 
061 
.046 
.046 
.031 
.031 
.016 
.016 
0056 
0054 
.0015 
.0014 
.00015 
00013 



rj range 



da /drj [nb] 



da/dr) [nb] 
forp* (D*) > 2.0 GeV 



-1.80 


-1.56 


1.19 


±0.14 


±8 


-1.56 


-1.32 


1.362 


±0.097 


±8 


-1.32 


-1.08 


1.418 


±0.071 


±° 

^0 


-1.08 


-0.84 


1.635 


±0.071 


±0 


-0.84 


-0.60 


1.629 


±0.069 


±° 

=t 


-0.60 


-0.36 


1.829 


±0.073 


=t 


-0.36 


-0.12 


1.731 


±0.071 


±° 


-0.12 


0.12 


1.878 


±0.081 


±° 
^o 


0.12 


0.36 


1.763 


±0.078 


=t 


0.36 


0.60 


1.927 


±0.090 


=t 


0.60 


0.84 


1.880 


±0.095 


=t 


0.84 


1.08 


2.025 


±0.097 


=t 


1.08 


1.32 


2.19 


±0.12 


±° 

=t 


1.32 


1.56 


1.97 


±0.17 


±° 
^o 


1.56 


1.80 


1.93 


±0.24 


=t 



0.460 


±0.078 


±1! 


0.500 


±0.051 


±8 


0.592 


±0.037 


it 


0.672 


±0.036 


± o 


0.728 


±0.038 


± () 


0.814 


±0.041 


±° 
^o 


0.836 


±0.042 


±° 
^o 


0.894 


±0.048 


^o 


0.824 


±0.044 


=t 


0.947 


±0.048 


± 


0.931 


±0.050 


±° 
^o 


0.939 


±0.049 


=t 


0.856 


±0.056 


=t 


0.764 


±0.077 


± 


0.876 


±0.107 


+° 



061 
047 
.050 
042 
.068 
043 
062 
.047 
.056 
.053 
.083 
.058 
.077 
064 
.070 
.067 
.065 
.060 
074 
.068 
.075 
.066 
.067 
.065 
.062 
.059 
.055 
.055 
.075 
(I'i'l 



p;J, range [GeV] 



dcr/dpj [nb/GeV] 



0.300 


0.700 


1.26 


±0.16 


± o 


0.700 


1.125 


1.83 


±0.14 


4-° 

±0 


1.125 


1.500 


2.22 


±0.15 


4-° 
±0 


1.500 


1.880 


2.39 


±0.14 


4-° 
=t 


1.880 


2.280 


2.02 


±0.11 


4-° 
=t 


2.280 


2.680 


1.417 


±0.086 




2.680 


3.080 


1.055 


±0.063 


4-° 
±0 


3.080 


3.500 


0.711 


±0.045 


4-° 


3.500 


4.250 


0.453 


±0.022 


4-° 


4.250 


6.000 


0.2028 


±0.0080 


4-° 
=t 


6.000 


11.000 


0.0287 


±0.0017 


±8 


11.000 


20.000 


0.00278 


±0.00062 


±8 



0.18 
0.18 



0.21 
0.18 
0.19 



z range 



da/dz [nb] 



5.12 


±0.59 


4_0.45 
^0.40 


9.42 


±0.59 


4_0.86 
^0.80 


10.36 


±0.48 


4_0.86 
^0.78 


9.66 


±0.41 


4_0.77 
^0.72 


9.30 


±0.36 


i 0.68 
^0.71 


4.97 


±0.16 


i_0.46 
^0.61 


1.086 


±0.082 


4_0.305 
^0.266 



dcr/dz [nb] 
for p* T {D*) > 2.0 GeV 



0.000 
0.100 
0.200 
0.325 
0.450 
0.575 
0.800 



0.100 
0.200 
0.325 
0.450 
0.575 
0.800 
1.000 



1.72 


±0.26 


4_0.59 

^o.eo 


4.52 


±0.29 


4_0.41 

^o^o 


5.44 


±0.24 


4_0.40 
^O^O 


4.91 


±0.21 


4_0.48 
=c 0.47 


3.71 


±0.16 


4.0.42 
=c 0.42 


1.156 


±0.066 


4.0.167 

^o.ies 


0.347 


±0.038 


4.0.077 

^o.oso 



Table 6: Differential D* cross section as a function of p T (D*), rj(D*), Pt(D*) and z(D*) in the 
kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7, \rj(D )\ < 1.8, p T (D*) > 1.25 GeV. 
The first quoted uncertainty is statistical and the second is systematic. 
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•q range 



d 2 a/d V dp T [nb/GcV] 



d 2 <r/dr)dp* T [nb/GeV] 



1.25 <p T < 2.00 GeV 



0.30 < p T < 1.25 GeV 



-1.8 
-1.2 
-0.6 
0.0 
0.6 
1.2 



-1.2 
-0.6 
0.0 
0.6 
1.2 
1.8 



0.760 
0.701 
0.704 
0.663 
0.760 
1.23 



±0.074 
±0.055 
±0.057 
±0.062 
±0.078 
±0.15 



±° 
± () 



058 
055 
052 
052 
052 
052 
049 
049 
053 
054 
09 



0.547 
0.453 
0.443 
0.398 
0.427 
0.56 



±0.051 
±0.035 
±0.038 
±0.046 
±0.084 
±0.18 



_i_0. 

=t o. 

iO. 

^o. 

iO. 
=t 0. 
iO. 

^o. 

iO. 

iO. 
± 0. 



041 
041 
037 
037 
051 
051 
046 
047 
072 
070 
04 

05 



2.00 < p T < 2.75 GeV 



1.25 < p T < 2.00 GeV 



-1.8 
-1.2 
-0.6 
0.0 
0.6 
1.2 



-1.2 
-0.6 
0.0 
0.6 
1.2 
1.8 



0.573 
0.586 
0.688 
0.703 
0.783 
0.723 



±0.042 
±0.030 
±0.034 
±0.043 
±0.045 
±0.062 



±° 

±° 
^o 

±° 

± 



044 
044 
042 
041 
051 
051 
050 
050 
055 
055 
051 
■ - 



0.481 

0.583 
0.640 
0.614 
0.708 
0.77 



±0.047 
±0.040 
±0.044 
±0.060 
±0.079 
±0.16 



^o. 
,o. 
=t o. 

4-0. 
^O. 
_l_0. 
31 0. 
_l_0. 

^o. 

iO. 
± 0. 



038 
035 
041 
041 
047 
045 
043 
048 
051 
057 
06 

06 



2.75 < p T < 4.00 GeV 



2.00 < p* T < 3.00 GeV 



-1.8 
-1.2 
-0.6 
0.0 
0.6 
1.2 



-1.2 
-0.6 
0.0 
0.6 
1.2 
1.8 



0.227 


±0.017 


1O.OI8 
^O.Ol? 


0.336 


±0.014 


_i_0.023 
=c 0.023 


0.359 


±0.014 


_i_0.025 
=c 0.025 


0.401 


±0.016 


+0.028 
=c 0.028 


0.377 


±0.017 


+0.027 
=c 0.026 


0.304 


±0.024 


+0.024 
3= 0.024 



0.336 
0.390 
0.392 
0.474 
0.549 
0.530 



±0.028 
±0.021 
±0.024 
±0.028 
±0.033 
±0.057 



,0. 
=t o. 



± o 

31 
31 
3= 



02!) 
029 
029 
027 
028 
028 
034 
035 
040 
040 
038 
039 



4.00 < p T < 6.00 GeV 



3.00 < p* T < 6.00 GeV 



1.8 


-1.2 


0.0368 


±0.0051 


± 


1.2 


-0.6 


0.1017 


±0.0051 


± 


0.6 


0.0 


0.1480 


±0.0059 


± 


0.0 


0.6 


0.1502 


±0.0068 


± 


0.6 


1.2 


0.1503 


±0.0068 


± 


1.2 


1.8 


0.0991 


±0.0095 


± 



0034 
0034 
0072 
0074 
0104 
0104 
0107 
0108 
0106 
0109 
0075 
0079 



0.0419 


±0.0047 


± 


0.0875 


±0.0042 


± 


0.1179 


±0.0048 


± 


0.1203 


±0.0055 


± 


0.1185 


±0.0063 


± 


0.0947 


±0.0095 


± 



0046 
0040 
0085 
0068 
0099 
0095 
0106 
0103 
0103 
0095 
0077 
0072 



6.00 <p T < 20.00 GeV 



6.00 < p T < 20.00 GeV 



1.8 


-1.2 


0.00073 


±0.00030 


± 


1.2 


-0.6 


0.00243 


±0.00035 


± 


0.6 


0.0 


0.00653 


±0.00050 


± 


0.0 


0.6 


0.00761 


±0.00053 


± 


0.6 


1.2 


0.00724 


±0.00057 


± 


1.2 


1.8 


0.00462 


±0.00076 


± 



00009 
00012 
00018 
00018 
00047 
00050 
00056 
00057 
00058 
00064 
00082 
000*2 



0.00030 


±0.00027 


± 


0.00196 


±0.00028 


± 


0.00377 


±0.00040 


± 


0.00439 


±0.00047 


± 


0.00571 


±0.00054 


± 


0.00257 


±0.00079 


± 



00004 
00005 
00021 
00016 
00036 
00029 
00051 
00044 
00052 
00051 
00020 
00021 



Table 7: Double differential D* cross sections as a function of rj(D*) and p T (D*) and as a 
function of r)(D*) and p* T {D*) in the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7, 
\rj(D*)\ < 1.8, p T (D*) > 1.25 GeV. The first quoted uncertainty is statistical and the second 
is systematic. 
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y range 



d 2 a/dydQ 2 [nb/GcV 2 ] 
forp T (D*) > 1.25 GeV, \t)(D*)\ < 



1.8 



d 2 a/dydQ 2 [nb/GeV 2 ] 
for p T (D*) > 1.5 GeV, \r](D*) 



< 1.5 



5 < Q 2 < 9 GeV 2 



0.02 


0.05 


2.27 


±0.19 


4_0.22 
^0.24 


0.05 


0.09 


1.95 


±0.16 


4_0.14 
^0.14 


0.09 


0.16 


1.767 


±0.096 


4_0.127 
^0.127 


0.16 


0.32 


0.954 


±0.052 


4_0.077 
^0.072 


0.32 


0.70 


0.361 


±0.024 


4_0.030 
^0.027 



1.23 


±0.14 


_i_0.11 

^o.n 


1.57 


±0.12 


_i_0.12 
=c 0.14 


1.378 


±0.077 


_i_0.114 

^o.m 


0.839 


±0.039 


_l_0.077 
^O.OSO 


0.243 


±0.018 


_l_0.024 
=c 0.025 



9 < - 



< 14 GeV 



0.02 
0.05 
0.09 
0.16 
0.32 



0.05 
0.09 
0.16 
0.32 
0.70 



0.845 
0.953 
0.687 
0.447 
0.193 



±0.101 
±0.085 
±0.054 
±0.029 
±0.015 



4_0.078 
^0.086 
4_0.072 
^0.075 
4_0.050 
^0.051 
4_0.032 
^0.031 
4_0.015 
^0.015 



0.457 

0.745 
0.581 
0.414 
0.144 



±0.082 
±0.068 
±0.044 
±0.028 
±0.011 



_i_0. 

± o. 

_i_0. 

± () 

i 0. 

^o 

±8: 
±8: 



045 
042 
065 
074 
048 
052 
031 
035 
013 
Ql:i 



14 < Q 2 < 23 GeV 2 



0.02 


0.05 


0.444 


±0.055 


± 


0.05 


0.09 


0.434 


±0.040 


± 


0.09 


0.16 


0.356 


±0.028 


± 


0.16 


0.32 


0.249 


±0.015 


± 


0.32 


0.70 


0.0887 


±0.0078 


± 



0.052 
0.049 
0.030 
0.033 
0.030 
0.029 
0.018 
0.018 
0.0079 
0.0077 



0.249 


±0.032 


± 


0.359 


±0.035 


± 


0.303 


±0.021 


± 


0.208 


±0.012 


± 


0.0659 


±0.0055 


± 



022 

023 

030 

033 

030 

032 

016 

019 

0065 

- 



23 < Q 2 < 45 GeV 2 



0.02 


0.05 


0.105 


±0.016 


± 


0.05 


0.09 


0.160 


±0.016 


± 


0.09 


0.16 


0.125 


±0.011 


± 


0.16 


0.32 


0.0885 


±0.0059 


± 


0.32 


0.70 


0.0375 


±0.0031 


± 



0.012 
0.011 
0.013 
0.014 
0.010 
0.009 
0.0062 
0.0064 
0.0031 
0.0030 



0.087 


±0.012 


± 


0.120 


±0.013 


± 


0.1211 


±0.0095 


± 


0.0744 


±0.0046 


± 


0.0304 


±0.0024 


± 



014 

014 

011 

012 

0107 

0120 

0056 

0065 

0027 



45 < i 



100 GeV 2 



0.02 


0.05 


0.0150 


±0.0085 


± 


0.05 


0.09 


0.0302 


±0.0054 


± 


0.09 


0.16 


0.0258 


±0.0034 


± 


0.16 


0.32 


0.0235 


±0.0022 


± 


0.32 


0.70 


0.0097 


±0.0011 


± 



0.0015 
0.0020 
0.0024 
0.0024 
0.0021 
0.0020 
0.0019 
0.0018 
0.0008 
0.0008 



0.0054 


±0.0026 


± 


0.0249 


±0.0041 


± 


0.0215 


±0.0028 


± 


0.0236 


±0.0019 


± 


0.00729 


±0.00085 


± 



0007 

0005 

0022 

0025 

0023 

0023 

0022 

0024 

00065 

00072 



Table 8: Double differential D* cross sections as a function of y and Q 2 in two differ- 
ent kinematic ranges: \r)(D*)\ < 1.8 and p T (D*) > 1.25 GeV or \r)(D*)\ < 1.5 and 
p T (D*) > 1.5 GeV. The first quoted uncertainty is statistical and the second is systematic. 



23 



Q 2 range [GeV 2 



da/dQ 2 [nb/GeV 2 



5.0 


6.0 


6.0 


8.0 


8.0 


10.0 


10.0 


13.0 


13.0 


19.0 


19.0 


27.5 


27.5 


40.0 


40.0 


60.0 


60.0 


100.0 



0.552 

0.398 

0.278 

0.1983 

0.1236 

0.0679 

0.0374 

0.01562 

0.00724 



±0.032 

±0.016 

±0.013 

±0.0088 

±0.0042 

±0.0028 

±0.0017 

±0.00095 

±0.00045 



_l_0.046 

=c 0.043 

_i_0.031 

^O.OSO 

_i_0.020 

=c 0.020 

_i_0.0143 

=c 0.0143 

_i_0.0088 

^o.ooss 

_i_0.0048 
=c 0.0048 
_i_0.0027 
=c 0.0027 
_i_0. 00113 
^O. 00110 
i 0.00053 
^O. 00051 



y range 



da/dy [nb] 



0.02 
0.05 
0.09 
0.13 
0.18 
0.26 
0.36 
0.50 



0.05 
0.09 
0.13 
0.18 
0.26 
0.36 
0.50 
0.70 



12.21 
16.39 
15.89 
12.71 
10.90 
6.85 
4.24 
2.13 



±0.64 
±0.69 
±0.72 
±0.61 
±0.42 
±0.31 
±0.22 
±0.17 



_i_1.47 
^1.47 
_i_1.21 
^l^O 
_i_1.30 
^1.25 
_i_0.92 
=c 0.92 
_i_0.84 
=c 0.82 
_i_0.53 

^o.si 

_i_0.37 
=c 0.34 
_i_0.20 
^Q.18 



x range 



da/dx [nb] 



0.00007 
0.00020 
0.00035 
0.00060 
0.00100 
0.00170 
0.00330 



0.00020 
0.00035 
0.00060 
0.00100 
0.00170 
0.00330 
0.05000 



3320 
4780 
3430 
2034 
1225 
446 
10.15 



±200 

±220 

±130 

±81 

±46 

±18 

±0.44 



_i_270 
=c 250 
_i_370 

^so 

_i_250 
=c 250 
_l_147 
=c 145 
_i_90 
^89 
,34 
± 34 
_i_0.84 
^Q.82 



Table 9: Differential D* cross section as a function of Q 2 , y and x in the kinematic range of 
5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7, \rj(D*)\ < 1.5 and p T (D*) > 1.5 GeV. The first quoted 
uncertainty is statistical and the second is systematic. 
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Pt range [GeV] 



da/dp T [nb/GeV] 



1.50 


1.88 


2 


34 


±0.15 


i 0.17 
^0.17 


1.88 


2.28 


2 


042 


±0.093 


i 0.144 
^0.144 


2.28 


2.68 


1 


959 


±0.070 


+0.140 
^O.MO 


2.68 


3.08 


1 


384 


±0.050 


+0.096 
^0.096 


3.08 


3.50 


1 


152 


±0.043 


+0.080 
^0.079 


3.50 


4.00 





814 


±0.028 


+0.056 
^0.056 


4.00 


4.75 





575 


±0.018 


+0.040 
^0.040 


4.75 


6.00 





2714 


±0.0088 


+0.0189 
^0.0187 


6.00 


8.00 





0851 


±0.0037 


+0.0058 
^0.0058 


8.00 


11.00 





0211 


±0.0015 


+0.0017 
^0.0016 


11.00 


20.00 





00178 


±0.00028 


+0.00013 
^O. 00012 


r\ r 


angc 


dtr/drj [nb] 


-1.50 


-1.25 


1 


229 


±0.077 


+0.090 
^0.088 


-1.25 


-1.00 


1 


319 


±0.062 


+0.098 
^0.094 


-1.00 


-0.75 


1 


501 


±0.061 


+0.113 
^0.108 


-0.75 


-0.50 


1 


635 


±0.065 


+0.118 
^0.116 


-0.50 


-0.25 


1 


569 


±0.063 


+0.112 
^0.109 


-0.25 


0.00 


1 


629 


±0.066 


+0.118 
^0.116 


0.00 


0.25 


1 


667 


±0.070 


+0.117 
^0.117 


0.25 


0.50 


1 


677 


±0.074 


+0.121 
^0.119 


0.50 


0.75 


1 


756 


±0.078 


+0.126 
^0.124 


0.75 


1.00 


1 


746 


±0.080 


+0.131 
^0.128 


1.00 


1.25 


2 


024 


±0.095 


+0.150 
^0.146 


1.25 


1.50 


1 


73 


±0.12 


+0.13 
=c 0.12 


range [GeV] 


da/dp* T [nb/GeV] 


0.300 


0.700 





75 


±0.13 


+0.13 
^0.13 


0.700 


1.125 


1 


34 


±0.12 


+0.14 
=l= 0.14 


1.125 


1.500 


1 


48 


±0.13 


+0.16 
^0.16 


1.500 


1.880 


1 


62 


±0.12 


+0.11 
=c 0.12 


1.880 


2.280 


1 


511 


±0.093 


+0.106 
^0.108 


2.280 


2.680 


1 


163 


±0.073 


+0.087 
^0.086 


2.680 


3.080 





884 


±0.055 


+0.061 
^0.061 


3.080 


3.500 





570 


±0.039 


+0.041 
^0.040 


3.500 


4.250 





403 


±0.020 


+0.030 
^0.030 


4.250 


6.000 





1785 


±0.0069 


+0.0166 
^0.0156 


6.000 


11.000 





0269 


±0.0015 


+0.0026 
^0.0025 


11.000 


20.000 





00186 


±0.00041 


+0.00020 
^O. 00016 



z range 



da/dz [nb] 



0.000 
0.100 
0.200 
0.325 
0.450 
0.575 
0.800 



0.100 
0.200 
0.325 
0.450 
0.575 
0.800 
1.000 



3.29 
7.02 
8.22 
7.59 
7.40 
4.06 
0.861 



±0.41 
±0.44 
±0.36 
±0.31 
±0.28 
±0.13 
±0.064 



+0.29 

^0.29 

+0.58 

^0.58 

+0.70 

^0.70 

+0.65 

^0.63 

+0.59 

^0.59 

+0.34 

^0.35 

+0.233 

^0.204 



Table 10: Differential D* cross section as a function of p T (D*), rj(D*), p* T (D*) and z(D*) in 
the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7, \rj(D*)\ < 1.5 and p T (D*) > 
1.5 GeV. The first quoted uncertainty is statistical and the second is systematic. 
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Q 2 [GeV 2 ] 


X 


HVQDIS 

Ff S ext [%] 


Sstat [%] 


Ssyst [%] 


CASCADE 

Ff S ext [%] 


6.5 


1.3 


10~ 4 


0.2160 


+8.5 
8.7 


±6 


7 


+ 7.7 
8 1 


0.2005 


+ 16.6 
7.3 


6.5 


3.2 


10~ 4 


0.1576 


+4.3 
^3.2 


±5 


5 


+ 7.7 
8 


0.1634 


+ 12.3 
^10.7 


6.5 
6.5 


5.0 
8.0 


10~ 4 

io- 4 


0.1516 
0.1036 


+4.2 

i 5.7 
±3.4 


±5 
±8 


4 
1 


±73 
+7.2 
±7.2 


0.1597 
0.1153 


+ 11.1 

^n.g 
+6.0 

=c 12.6 


6.5 


2.0 


IO -3 


0.0735 


+ 10.8 
^7.2 


±8 


6 


+9.9 
^lO^ 


0.1044 


+7.2 
±12.2 


12.0 


3.2 


10~ 4 


0.2829 


+8.7 
±5.6 


±7 


7 


+ 7.6 
7 9 


0.2727 


+ 15.5 
^8.6 


12.0 


5.0 


10~ 4 


0.2123 


,3.1 
2 9 


±6 


6 


+ 7.1 
7 1 





2169 


+ 12.5 
10 7 


12.0 


8.0 


10~ 4 





1689 


_l_4.6 
±2.3 


±7 


8 


+7:4 





1779 


+io!o 


12.0 


2.0 


IO" 3 





1226 


.6.1 
±3.5 


±8 


9 


+ 7.7 
7.7 





1353 


+8.7 
=c 14.4 


12.0 


3.2 


IO" 3 





0773 


+ 11.6 
±7.4 


±12 





+9.6 
^9.8 





1125 


+5.7 


18.0 


5.0 


10~ 4 





3221 


+4.6 
±5.0 


±8 


8 


+8.5 
±9.1 





3045 


+ 15.2 
±9.2 


18.0 


8.0 


io~ 4 





2899 


+3.8 
±2.1 


±6 


1 


+7.3 
±7.2 





2964 


+ 11.7 

^lo.e 


18.0 
18.0 


1.3 
3.2 


IO" 3 
IO" 3 






2167 
1368 


+4.0 
=c 2.9 
+5.3 
±3.5 


±8 
±9 




3 


+8.2 
^8.3 
+7.4 
±7.2 






2202 
1471 


+11.0 

=c 12.9 
+ 10.1 

^le.o 


18.0 


5.0 


IO" 3 





1033 


+ 13.6 
^6.0 


±12 


5 


+11.3 

^n.s 





1455 


+8.7 
^ls.e 


35.0 
35.0 


8.0 
1.3 


io~ 4 

IO" 3 






3958 
3188 


+3.6 
^3.0 
+2.8 
^2.4 


±8 
±6 


3 
7 


+8.0 
^8.2 
+ 7.2 
±7.1 






3620 
3092 


+ 14.0 
^11.7 
+ 11.9 
±13.5 
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Table 11: F 2 CC in bins of Q 2 and x extracted from measured D* cross sections with two dif- 
ferent programs, HVQDIS and CASCADE. The extrapolation uncertainty 5 ext is determined 
by varying model parameters within a program. The statistical (S sta t) and systematic (8 sys t) 
uncertainties arise from the determination of the D* cross section and are the same for both 
programs. 
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Figure 1: Distribution of the mass difference Am = m(fC T 7r ± 7r^ = ) — m(_ft' =F 7r ± ) for D* can- 
didates in the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7, \r](D*)\ < 1.8 and 
Pt(D*) > 1.25 GeV. The histogram shows the wrong charge combinations, K ± tt ± tcT. The 
solid line represents the result of the fit described in the text. 
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Figure 2: Differential D* cross section as a function of the photon virtuality Q 2 , the inelasticity 
y and Bjorken x. The measurements correspond to the kinematic range of 5 < Q 2 < 100 GeV 2 , 
0.02 < y < 0.7, \rj(D*)\ < 1.8 and p T (D*) > 1.25 GeV. The data are shown as points, 
the inner error bars show the statistical error, the outer error bars represent the statistical and 
systematic errors added in quadrature. The data are compared to predictions by the MC program 
RAPGAP with two different proton parton densities and by the MC program CASCADE. In 
the lower part of the figures the normalised ratio _R norm of theory to data (equation 0) is shown, 
which has reduced normalisation uncertainties. 
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Figure 3: Differential D* cross section as a function of the photon virtuality Q 2 , the inelasticity 
y and Bjorken x. The measurements correspond to the kinematic range of 5 < Q 2 < 100 GeV 2 , 
0.02 < y < 0.7, \tj(D*)\ < 1.8 and p T {D*) > 1.25 GeV. The data are shown as points, the 
inner error bars show the statistical error, the outer error bars represent the statistical and sys- 
tematic errors added in quadrature. The data are compared to predictions by the next-to-leading 
order calculation HVQDIS with two different proton parton densities. The bands indicate the 
theoretical uncertainties (table |2]). In the lower part of the figures the normalised ratio _R norm of 
theory to data (equation |3]) is shown, which has reduced normalisation uncertainties. 
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Figure 4: Differential D* cross section as a function of the transverse momentum p T (D*) 
and pseudo-rapidity r](D*) in the laboratory frame, the transverse momentum p^(D*) in the 
7P centre-of-mass frame and the D* inelasticity z(D*). The measurements correspond to 
the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7 and \rj(D*)\ < 1.8 and 
Pt(D*) > 1.25 GeV. The data are shown as points, the inner error bars show the statistical 
error, the outer error bars represent the statistical and sy stematic errors added in quadrature. 
The data are compared to predictions by the MC program RAPGAP with two different proton 
parton densities and by the MC program CASCADE. In the lower part of the figures the nor- 
malised ratio R norm of theory to data (equation |3]) is shown, which has reduced normalisation 
uncertainties. 
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Figure 5: Differential D* cross section as a function of the transverse momentum pt(D*) 
and pseudo-rapidity r)(D*) in the laboratory frame, the transverse momentum p^(D*) in the 
7P centre-of-mass frame and the D* inelasticity z(D*). The measurements correspond to 
the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7 and \r)(D*)\ < 1.8 and 
Pt(D*) > 1.25 GeV. The data are shown as points, the inner error bars show the statistical 
error, the outer error bars represent the statistical and sy stematic errors added in quadrature. 
The data are compared to predictions by the next-to-leading order calculation HVQDIS with 
two different proton parton densities. The bands indicate the theoretical uncertainties (table |2]). 
In the lower part of the figures the normalised ratio _R norm of theory to data (equation |3]) is 
shown, which has reduced normalisation uncertainties. 
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Figure 6: Double differential D* cross section as a function of the transverse momentum 
Pt(D*) and pseudo-rapidity r](D*) in the laboratory frame. The measurements correspond 
to the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7 and \rj(D*)\ < 1.8 and 
Pt(D*) > 1.25 GeV. The data are shown as points, the inner error bars show the statistical 
error, the outer error bars represent the statistical and sy stematic errors added in quadrature. 
The data are compared to predictions by the MC program RAPGAP with two different proton 
parton densities and by the MC program CASCADE. 
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Figure 7: Double differential D* cross section as a function of the transverse momentum 
Pt(D*) and pseudo-rapidity r](D*) in the laboratory frame. The measurements correspond 
to the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7 and \rj(D*)\ < 1.8 and 
Pt(D*) > 1.25 GeV. The data are shown as points, the inner error bars show the statistical 
error, the outer error bars represent the statistical and sy stematic errors added in quadrature. 
The data are compared to predictions by the next-to-leading order calculation HVQDIS with 
two different proton parton densities. The bands indicate the theoretical uncertainties (table H). 
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Figure 8: Double differential D* cross section as a function of the transverse momentum in the 
7P centre-of-mass frame p^(D*) and the pseudo-rapidity r](D*) in the laboratory frame. The 
measurements correspond to the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7 
and \r)(D*)\ < 1.8 and pt(D*) > 1.25 GeV. The data are shown as points, the inner error 
bars show the statistical error, the outer error bars repres ent the statistical and systematic errors 
added in quadrature. The data are compared to predictions by the MC program RAPGAP with 
two different proton parton densities and by the MC program CASCADE. 
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Figure 9: Double differential D* cross section as a function of the transverse momentum in the 
7P centre-of-mass frame p^(D*) and the pseudo-rapidity i](D*) in the laboratory frame. The 
measurements correspond to the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7 
and \i](D*)\ < 1.8 and pr(D*) > 1.25 GeV. The data are shown as points, the inner error 
bars show the statistical error, the outer error bars repres ent the statistical and systematic errors 
added in quadrature. The data are compared to predictions by the next-to-leading order calcu- 
lation HVQDIS with two different proton parton densities. The bands indicate the theoretical 
uncertainties (table |2]). 
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Figure 10: Double differential D* cross section as a function of the photon virtuality Q 2 and the 
inelasticity y. The measurements correspond to the kinematic range of 5 < Q 2 < 100 GeV 2 , 
0.02 < y < 0.7 and \r](D*)\ < 1-8 and p T {D*) > 1.25 GeV. The data are shown as points, 
the inner error bars show the statistical error, the outer error bars represent the statistical and 
systematic errors added in quadrature. The data are compared to predictions by the MC program 
RAPGAP with two different proton parton densities and by the MC program CASCADE. 
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Figure 11: Double differential D* cross section as a function of photon virtuality Q 2 and the 
inelasticity y. The measurements correspond to the kinematic range of 5 < Q 2 < 100 GeV 2 , 
0.02 < y < 0.7 and \ij(D*)\ < 1.8 and p T {D*) > 1.25 GeV. The data are shown as points, the 
inner error bars show the statistical error, the outer error bars represent the statistical and sys- 
tematic errors added in quadrature. The data are compared to predictions by the next-to-leading 
order calculation HVQDIS with two different proton parton densities. The bands indicate the 
theoretical uncertainties (table 0). 
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Figure 12: Differential D* cross section as a function of photon virtuality Q 2 , the inelasticity y 
and Bjorken x. The measurements correspond to the kinematic range of 5 < Q 2 < 100 GeV 2 , 
0.02 < y < 0.7, \r](D*)\ < 1-8, Pt{D*) > 1.25 GeV with an additional cut on the D* trans- 
verse momentum in the "fp centre-of-mass frame p^(D*) > 2.0 GeV. The data are shown as 
points, the inner error bars show the statistical error, the outer error bars represent the statistical 
and systematic errors added in quadrature. The data are compared to a prediction to next-to- 
leading order in the ZM-VFNS and to HVQDIS. The bands indicate the theoretical uncertainties 
(table©. 
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Figure 13: Differential L>* cross section as a function of the transverse momentum p T (D*) and 
pseudo-rapidity T)(D*) in the laboratory frame and the D* inelasticity z(D*). The measurements 
correspond to the kinematic range of 5 < Q 2 < 100 GeV 2 , 0.02 < y < 0.7, \r){D*)\ < 1.8, 
Pt(D*) > 1.25 GeV with an additional cut on the D* transverse momentum in the centre- 
of-mass frame p^(D*) > 2.0 GeV. The data are shown as points, the inner error bars show 
the statistical error, the outer error bars represent the statistical and systematic errors added in 
quadrature. The data are compared to a prediction to next-to -leading order in the ZM-VFNS 
and to HVQDIS. The bands indicate the theoretical uncertainties (tabled. 
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D* in DIS 




Figure 14: Differential D* cross section as a function of the photon virtuality Q 2 . The 
measurements correspond to the kinematic range of 0.02 < y < 0.7, \r](D*)\ < 1.5 and 
Pt(D*) > 1.5 GeV. The data of this measurement (points) are shown in a phase space with 
stronger restrictions on r](D*) and pt(D*) to be comparable to a previous measurement at 
higher Q 2 lfT5l (triangles). The inner error bars show the statistical error, the outer error bars 
represent the statistical and systematic errors added in quadrature. The data are compared to 
predictions by the MC program RAPGAP with two different proton PDFs and by the MC pro- 
gram CASCADE (left) and to predictions by the next-to-leading order calculation HVQDIS 
with two different proton PDFs (right). 
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Figure 15: Extrapolation factors from the visible phase space (table© to the total phase space 
for the D* meson as determined from HVQDIS and CASCADE. The error bars show the ex- 
trapolation uncertainty which is determined by varying the theory parameters listed in table |2] 
for HVQDIS and in tabled] for CASCADE. 
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Figure 16: F 2 CC as derived from D* data with HVQDIS (points). The inner error bars show the 
statistical uncertainty, the outer error bar the statistical and experimental systematic uncertainty 
added in quadrature. The extrapolation uncertainty within the HVQDIS model is shown as 
blue band in the bottom of the plots. The outer (orange) band shows the model uncertainty 
obtained from the difference in F 2 CC determined with HVQDIS and CASCADE. The data are 
compared to the measurement of F 2 CC with the HI vertex detector [14] (open squares), to NLO 
DGLAP predictions from HVQDIS with two different proton PDFs, and to the F 2 CC prediction 
of HERAPDF1.0. 
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Figure 17: F 2 CC as derived from D* data with HVQDIS (points). The inner error bars show the 
statistical uncertainty, the outer error bar the statistical and experimental systematic uncertainty 
added in quadrature. The extrapolation uncertainty within the HVQDIS model is shown as blue 
band in the bottom of the plots. The outer (orange) band shows the model uncertainty obtained 
from the difference in determined with HVQDIS and CASCADE. The data are compared 
to the measurement of F| c with the HI vertex detector ffT4l (open squares) and to predictions 
from the global PDF fits CT10 (dashed line), MSTW08 at NNLO (dark dashed-dotted line), 
NNPDF2.1 (shaded band) and ABKM (light dashed-dotted line). 
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Figure 18: F 2 CC as a function of Q 2 for different x, as derived from D* data with HVQDIS 
(points). The inner error bars show the statistical uncertainty, the outer error bar the total un- 
certainty, including statistical, experimental systematic, extrapolation and model uncertainty 
added in quadrature. The data are compared to the measurement of F 2 CC with the HI vertex 
detector [14] (open squares), to NLO DGLAP predictions from HVQDIS with two different 
proton PDFs, and to the Ff prediction of HERAPDF1 .0. 
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